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Abstract. We present gauge invariant spectral Cauchy characteristic extraction. 
We compare gravitational waveforms extracted from a head-on black hole merger 
simulated in two different gauges by two different codes. We show rapid convergence, 
demonstrating both gauge invariance of the extraction algorithm and consistency 
between the legacy Pitt null code and the much faster Spectral Einstein Code 
(SpEC). 



1. What is CCE? What is gravitational waveform gauge invariance? 



Figure 1: Cauchy characteristic extraction. A Cauchy evolution of the Einstein field 
equation proceeds on a space-like foliation (green). A finite spheroidal worldtube E at 
areal radius forms the inner boundary to a characteristic evolution on a null foliation 
(red). Based on a coordinate system derived from the Cauchy data on the worldtube, 
gravitational information is propagated to compactified future null infinity . At an 
inertial coordinate system is co-evolved, in which the desired gauge-invariant waveform 
can be expressed. 

The strong gravitational radiation produced in the inspiral and merger of binary 
black holes has been a dominant motivation for the construction of gravitational 
wave observatories. The details of the gravitational waveform supplied by numerical 
simulation is a key theoretical tool to fully complement the sensitivity of the LIGO, 
Virgo, GEO, and KAGRA observatories, by enhancing detection and providing useful 
scientihc interpretation of the gravitational signal[Il El |3l ll]. Gharacteristic evolution 
coupled to Gauchy evolution via Gauchy-characteristic extraction (GGE) provides the 
most accurate numerical computation of the Bondi news function, which determines 
both the waveform and the radiated energy and momentum at null inhnity. 








In CCE, the Cauchy evolution is used to supply boundary data on a timelike 
inner worldtube necessary to carry out a characteristic evolution extending to future 
null inhnity where the radiation is computed using the geometric methods 

developed by Bondi et a/. [5], Sachs [6], and Penrose |7], as depicted in Fig. 
More intuitive methods, including intrinsically inertial compactihed hyperboloidal 
formulations [HI IHl HD] have not yet found adoption in the evolution of binary black 
holes. CCE is an initial-boundary value problem based upon a timelike worldtube cu¬ 
lt has been implemented as a characteristic evolution code, the Pitt null code [121CH], 
which incorporates a Penrose compactihcation of the spacetime, and which has 
subsequently been extended to higher order methods by Reisswig et al. na. It has 
more recently been implemented as a spectral code within the Spectral Einstein Code 
(SpEC) by Handmer and Szilagvi[TH]. upon which the present work is based. 

One technical complication introduced by CCE is that the coordinates induced 
on by the computational Cauchy coordinates on the inner worldtube do not 
correspond to inertial observers, i.e., to the coordinates intrinsic to a distant freely 
falling and non-rotating observatory. The gravitational waveform Erst obtained 
in the “computational coordinates” of CCE is in a scrambled form. This gauge 
ambiguity in the waveform is removed by constructing the transformation between 
computational coordinates and inertial coordinates at . There still remains the 
freedom in the choice of inertial observers. In special relativistic theories, this freedom 
is reduced to the translations and Lorentz transformations of the Poincare group. As 
explained in Sec.[^ in an asymptotically flat space time the corresponding asymptotic 
symmetry group consists of supertranslations and Lorentz transformations. This 
freedom governs the redshift and initial phase of the waveform. 

A physically relevant calculation of the radiation flux must also be referred 


to such inertial coordinates at . In this paper, the calculation of the energy- 
momentum flux via the Bondi news function is hrst carried out in the induced 
worldtube coordinates and then transformed to the inertial coordinates. 

2. Characteristic Formalism 

The characteristic formalism is based upon a family of outgoing null cones emanating 
from an inner worldtube and extending to inhnity where they foliate into spherical 
slices. We let u label these hypersurfaces, {A = 2, 3) be angular coordinates that 
label the null rays, and r be a surface area coordinate along the outgoing null cones. 

Employing the conventions used in [15], in the resulting = {u,r,y^) 
coordinates, the metric takes the Bondi-Sachs form 

ds^ = - {e^^{rW + 1) - r^hAsU^U^) du^ 

— 2e^^dudr — 2r‘^hABU^dudy^ + r^hAsdy^dy^ , (2.1) 

where h^^hsc = det(hAB) = det{qAB), with qab a unit sphere metric. In 

analyzing the Einstein equations, we also use the intermediate variable 

Qa = r^e~‘^^hABU^ . ( 2 . 2 ) 

The metric coefficients IT, Hab, Qa, /3 represent respectively the mass aspect, 
the spherical 2-metric, the shift and its radial derivative, and the lapse. The vector 
and tensor helds hAB-,U^,QA are expressed as spin-weighted helds by contracting 
them with the complex dyad for the unit sphere metric satisfying q'^qA = 0, 
q^qA = 2, q^ = q^^qB, with q^^qBC = and qAB = {qAqB + qAqB)/‘2. Under 
this convention, the spin-weighted functions U = U^qA and Q = QAq^, while 
J = hABd^Q^/‘^ uniquely determines the spherical 2-metric component of the general 
4-metric p. We chose a dyad consistent with the computational formulation of the 


spin-weight raising S operator [T^, given by = (—1, —i/ sin0) in standard spherical 
coordinates This is regular everywhere except the poles, which we can avoid 

through careful choice of grid points. It is worth noting that any choice of angular 
coordinates are possible. Other conventions use multiple patches to avoid singularities 
at the poles. 

A key feature of the Bondi-Sachs formulation is that the Einstein equations can 
be integrated along the outgoing characteristics in a sequential order. We use a form 
which hrst appeared in im and was implemented as the Pitt code in [T31 HH]: 


(2.3) 

{r^Q),r = - r^dJ + dK),r +2r%r-^(]),r +Nq , (2.4) 

U,r = r-^e^^Q + Nu , (2.5) 

{r^W),r = - 1 - + 8t/))„ +Nw , (2.6) 

and the evolution equation 


2(r J),„. = ((1 + rW)irJ),r ) -r-\r^m),r +2r-^e^Q^e^ - (rW)„ J + Nj , (2.7) 

where 

n = 2K- mK + ^(8V + 8^J) + ^(8^8^ - 8J8J) , (2.8) 

is the curvature scalar associated with hAs, = 1 + JJ and Nq, Nw, Nj are 
nonlinear terms given in [T3j. 

On each constant u hypersurface of the spacetime foliation, these equations are 
integrated in turn. Given J, /3 is solved, then U, Q, and W in turn, enabling the 
computation of J,„. J,„ permits a step forward in time and J is thus dehned on the 
next hypersurface. The radial compactihcation of inhnity is given by 

r = r„,p/(l-p) , i<p<l, 


(2.9) 


where the compactification parameter y^) is the (not necessarily constant) areal 

radius coordinate on the worldtube. 

Angular derivatives are implemented using the action of the 9 operator on spin- 
weighted spherical harmonics, e.g., Sf/ = q^q^UA.Bi where a colon denotes the 
covariant derivative with respect to qAB [IS]- In spherical coordinates, this takes 
the explicit form for a spin-weight-s field y 

Sv=-{An-e) (k + . (2,10) 

8o = -(sin-9) ^ A j (sin-e,) . (2,11) 

8 is the associated spin-weight lowering operator. 

The spectral algorithm used to solve these equations and the treatment of the 
nonlinear terms Nq, Nu, Nw, Nj are detailed in Handmer and Szilagyi[T5]. Here, 
we extend the characteristic spectral algorithm to calculating the gauge invariant 
Bondi news at 


3. Waveforms at 


For technical simplicity, the theoretical derivation of the waveform at inhnity is best 
presented in terms of an inverse surface-area coordinate i = 1/r, where £ = 0 at 
In the resulting = (u, i, conformal Bondi coordinates, the physical spacetime 
metric has the conformal compactihcation where is smooth at 


and, referring to the metric (2.1), takes the form|llj 


+ iW) - HabU^U^) du^ + 2e^^dudi (3.1) 

— 2hABU^dudx^ -I- hABdx^dx^ . 


As described in [191 [20|, both the Bondi news function N{u, x^) and the Newman- 





Penrose Weyl tensor component [21] 


= lim r'04 , (3.2) 

r^oo 

which describe the waveform, are determined by the asymptotic limit at of the 
tensor field 

. (3.3) 

This limit is constructed from the leading coefficients in an expansion of the metric 
about in powers of i. We thus write 

hAB = HAB+^CAB + 0[e) . (3.4) 

Conditions on the asymptotic expansion of the remaining components of the metric 
follow from the Einstein equations: 

/3 = if + 0(0, (3.5) 

U^ = L^ + 2ie^^H^^DBH + 0{f) , (3.6) 

and 

W = DaL^ + i{e‘^^n/2 + DaD^c^^ - 1) + , (3.7) 

where H and L are the asymptotic limits of f3 and U and where 77 and Da are the 
2-dimensional curvature scalar and covariant derivative associated with Hab- 
The expansion coefficients H, Hab, cab, and (all functions of u and 
completely determine the radiation field. One can further specialize the Bondi 
coordinates to be inertial at i.e., have Minkowski form, in which case H = 
= 0, Hab = Qab (the unit sphere metric) so that the radiation field is completely 
determined by cab- However, the characteristic extraction of the waveform is carried 
out in computational coordinates (determined by the Cauchy data on the extraction 
worldtube) so this inertial simplification cannot be assumed. 



In order to first compute the Bondi news function in the g^i, computational 
frame, it is necessary to determine the conformal factor u relating Hab to a unit 
sphere metric Qab, i-e., to an inertial conformal Bondi frame m satisfying 

Qab = ^‘^Hab ■ (3.8) 


(See [22] for a discussion of how the news in an arbitrary conformal frame is related 
to its expression in this inertial Bondi frame.) We can determine u by solving the 


elliptic equation governing the conformal transformation of the curvature scalar (2.8) 
to a unit sphere geometry: 


n^2(ui‘ + H-"^DaDb logw) . 


(3.9) 


The elliptic equation (3.9) need only be solved at the initial time where, with initial 
data J\j^+ = 0, DaDb simplihes to the 2-Laplacian on the unit sphere. Then, as 
described in the next section, application of the Einstein equations on determines 
the time dependence of u according to 


2h"9aloga; = —e DaL^ , 


(3.10) 


where h" = pi is the null vector tangent to the generators of We use (3.10) 


to evolve u along the generators of given a solution of (3.9) as initial condition. 


First recall some basic elements of Penrose compactihcation. In a general 
conformal frame with metric g^i, = VL^g^Bv, where = 0 on the vacuum Einstein 
equations = 0 take the form 

- g^,^ (2f2V^Vpf2 - 3(V^f2Vpf2) J = 0 . (3.11) 


It immediately follows that 


(V^fl)Vpf2U+ = 0 , 


(3.12) 






so that is a null hypersurface and that 


V^Vpff] U+ = 0 . 


(3.13) 


With respect to this frame, the construction of an inertial conformal frame 
proceeds as follows. We introduce a new conformal factor = caff, with = uj^g^v 


by requiring, in accord with (3.10), 


\2ffdaOJ + caVo-h'^] |jr+ = 0 , = g^'^WpVt 


(3.14) 


As a result, it follows from a straightforward calculation that 

V'’VpfiU+ = 0 , (3.15) 

i.e., in the g^^ conformal frame is null, shear- and divergence-free. 

It also follows that 

h‘"V<,h'^|,y+ = 0 , (3.16) 

where = gf’^'WpVL, i.e., in the g^^ frame, is an affinely parametrized null generator 
of ^+. 

To construct inertial coordinates {u,x^) on we hrst assign angular 

coordinates to each point of the initial spacelike spherical slice m = mq of 
We then propagate these coordinates along the generators of according to 

hPdpX^\j!+ = uj~^n'^dpX^\y+ = 0 . (3-17) 


In addition, we require 


n'^dpu\jr+ = uj ^h^dpu\j!+ = 1 , 


(3.18) 


so that u is an affine parameter along the generators in the g^v conformal frame. 



4. News 


The Bondi news function N is computed in the computational coordinates with 
the appropriate conformal transformation. It is then interpolated onto the inertial 
coordinates. The formalism follows that of [13], Appendix B (with a sign error in S 3 
corrected): 

N = ^—— ^Si + S 2 + “ (pU + S 3 — 4uj ^84 + 2cu ^S 5 j , (4-1) 



where A = and the Si terms are 


51 — + JJJ/u ~ 2 J KKgu ~ 2 J+ 2J^iu)/{K + 1 ) , 

5 2 = {QJ^eJJU + 2QJ^eKU + 2QJ^eU + QJ^eJ^U 

- 2dK,iJKU - 2dK^eJU + 2mjJK^i - 2WJJ^eK 

- 2WJJ,e + AdUKK^e + AWK^i + 2WJJJ,e - 2dUJKK^e 

- 2WJK^^ + AWJ^eK + AQUJ^e + QJ^eJJU + 2dJ^eKU + 25 
+ dJ^eJ'^U - 2dK^iJKU - 2QK^eJU + 2WJ‘^J^t - 2mjKK^i 

- 2mjK^i + 2Wj‘^K^i - 2WJJ^iK - 2mjJ,e)/{2{K + 1 )) , 

5 3 = -{J‘^Je + JJJ,e - 2JKK^e - 2JK^e + 2J^eK + 2J^e)/{K + 1 ) , 

5 4 = {QAQuJJ + 2 QAQujK + 25A5a; — QAQuJK 

- dABuJ - dudAJK - dudAJ + BABuJ^)/{2{K + 1)) , 

5 5 = {2d^AJJ + AQ'^AK + 45 M + 2Q‘^AJ‘^ - AddAJK 

- AddAJ + dAdJJP + 2dAdJJK + 25A5JJ + dAdJJ^J 

+ 2dAdJJK + 25A5 J J - 2QAQKJJK - AdAdKJJ - AS ASK K 

- ASASK - SASJJJK + 2SASJK + 25A5J - SAdJJ^K 
+ 2SASKJ^J - SJSAJJK - 2SJSAJJ - 2SJSAK 

- 2dJdA - SJSA.PK - 2SJSA.P + 2SKSA.PJ 
+ ASKSAJK + ASKSAJ + SASJJ^J 

+ SASJ.f - 2SASK.fK)/{A{K + 1)) . (4. 

In our implementation, derivatives are derived from spectrally calculated 
derivatives using the appropriate Jacobian. 



5. Results 


In our comparison tests of CCE, the worldtube boundary data were extracted from 
a simulation of an equal mass non-spinning head-on black hole collision, with initial 
separation of SOM. The control run (Isotropic) utilized the standard harmonic gauge 
damping identical to that in [23] throughout the head-on merger and ring-down. 
Harmonic gauge damping adds a dissipative forcing term to the wave equations 
satished by the harmonic Cartesian spatial coordinates {x, y, z). In order to diminish 
the effects of a custom designed gauge, we also compare with results of another run 
(HytZero) which turns off gauge damping in the harmonic ^/-direction, transverse to 
the x-direction motion of the black holes. These two high-resolution runs were used 
as boundary data for all the subsequent CCE runs. These runs include 3 different 
resolutions, 2 different codes, 2 different gauges, and 3 different extraction radii, for 
a total of 36 runs. 

As described in [15] , the SpEC characteristic evolution algorithm exploits spectral 
methods and innovative integral methods that greatly improve upon the speed and 
accuracy of the Pitt null code. This is seen as essential for for taking advantage of the 
efficiency of SpEC Cauchy evolution. The necessary improvement in efficiency has been 
preserved in the SpEC extraction module, as displayed in Table The comparison 
runs were performed using the current version of the Pitt code [2D] , which forms part 
of the Einstein Toolkit. 

The initial conditions and extraction parameters were deliberately chosen as a 
stressful test of the algorithms. In particular, at the beginning of the run the black 
hole excision boundary extends out to Cartesian radius R = 16M, which is very close 
to our smallest choice of extraction radius at i? = 30M. At this radius, gauge effects 
are highly signihcant and would make perturbative extraction schemes meaningless. 


in accordance with our intentions. One consequence of such an extreme choice is 
that differences between the Pitt and SpEC inertial frame and worldtube initialization 
procedures lead to noticeably different waveforms. Worldtube initialization involves 
supplying the “integration constants” from the Cauchy code, which allows radial 


integration of the characteristic hypersurface and evolution equations (2.3) - (2.7) 
from the worldtube to . In both Pitt null code and SpEC, the initial condition 
on J is determined by the inner boundary value, supplied by the Cauchy evolution, 
with a smooth roll off to zero at 

The extraction worldtube P is determined by a surface of constant Cartesian 
radius R. In the Pitt CCE code, the areal radius r^t of P lies between two surfaces 
of constant Cartesian radii Ri < r^t < R 2 and this carries over to the compactihed 
radial coordinate. As a result, interpolation is necessary to supply the integration 
constants, which introduces numerical error. In the SpEC CCE code, this interpolation 


error is avoided by introducing the compactihed radial coordinate (2.9), with range 
1/2 < p < 1 between P and 

Worldtube data from each run were extracted using both Pitt and SpEC CCE, at 
three different Cartesian radii: R = SOM, R = lOOM and R = 250M, as illustrated 
by the news function waveforms in Figs. and[^ respectively. In these hgures, the 
HytZero and Isotropic waveforms are so close that they appear on top of one another. 
The major discrepancy between the Pitt and SpEC waveforms is due to the worldtube 
interpolation error in the Pitt code. This is especially evident at small extraction 
radii, where there is strong “junk” radiation near the worldtube, which is inherent in 
the initial Cauchy data and its mismatch with the initial characteristic data. 


This interpolation error in the Pitt code converges away at larger radii, where 
the held gradients between Ri and R 2 become smaller. This is seen in Figs. and 





Run 

Pittl 

Pitt2 

Pitt3 

SpECl 

SpEC2 

SpEC3 

Nr 

100 

150 

200 

10 

12 

14 

N or L 

40 

60 

80 

12 

14 

17 

At/M 

0.1 

0.0666... 

0.05 

1.0 

0 . 666 ... 

0.5 

T (CPU hours) 

173 

274 

374 

0.7 

1.9 

3.1 


Table 1: Resolution parameters used for code convergence comparisons, with time steps 
At. Nr represents the radial grid sizes. The Pitt null code uses two stereographic patches 
with 2N‘^ total number of angular grid points. The SpEC code has total angular grid 
points. T is the CPU time taken for R = SOM, tfinal = 450M runs in the Isotropic gauge, 
and is representative for the other runs. All resolutions and codes were run from the same 
initial data. 

where the relative difference between the Pitt and SpEC news function waveforms 
is compared with the relative numerical error implied by convergence tests. 

Each run was computed at 3 different resolutions to monitor convergence, as 
indicated in Table In the following subsections, we hrst show convergence and the 
removal of gauge effects, separately for the Pitt and SpEC codes. Next, we compare 
T 4 waveforms and establish further agreement between the two codes. Finally, we 
examine the evolution of the inertial coordinates at relative to the worldtube 
coordinates induced by the Cauchy evolution. 

Comparison of the relative error Erei between dataset A and dataset B is 
computed according to 

Erel = logio 

where in convergence tests B is the highest resolution dataset, and the real parts of 
the {i, m) = ( 2 , 2 ) spherical harmonic modes are compared. 








Figure 2: Waveforms of the real part of the (2, 2) spherical harmonic mode of the news 
function, as computed by the Pitt null code and SpEC with extraction worldtube at 
R = SOM. Different initialization procedures at the worldtube give rise to a difference 
between the Pitt and SpEC waveforms, which is most pronounced at this small extraction 
radius. The different gauge choices. Isotropic and HytZero, do not have noticeable effect on 
this scale, indicating successful gauge effect removal in both codes. 


5.1. Pitt code convergence and removal of gauge effects 


Here, in order to establish a baseline, we examine the self convergence of the Pitt 
code for each of the extraction radii, using the three resolutions (Pittl,Pitt2,Pitt3) 
indicated in Table In Figs. | 


and 10, we see in both the Isotropic and HytZero 
gauges that the news function converges over the entire run. Indeed, Isotropic (solid 
lines) and HytZero (dashed lines) overlap completely. The hgures also plot the 
relative error in the news computed in both gauges, which is consistently below the 
numerical error implied by convergence tests for extraction worldtubes at R = SOM 
and R = lOOM. This verihes that the Pitt code successfully removes gauge effects. 
Furthermore, the hgures plot the relative error between the news computed in the 
















Figure 3: Graphs showing the relative difference between the real part of the (2,2) mode 
Pitt and SpEC news function waveforms for extraction radius R = SOM, in comparison to 
the relative numerical error implied by convergence tests, corresponding to the waveforms 
in Fig. While both SpEC and Pitt have comparable and consistent levels of error, the 
codes do not agree within that level of error at this extraction radius. 


worldtube coordinates and the inertial coordinates at . In the R = SOM case 
shown in Fig. the initial discrepancy is high due to the strong gauge effects of junk 
radiation. It does not fall below the relative error between the Isotropic and HytZero 
gauges until well after the signal has passed. This confirms that the transformation to 
inertial coordinates is essential for correctly removing gauge effects from the waveform. 
For extraction aX R = lOOM shown in Fig. the relative error between worldtube 
and inertial coordinates has dropped below the Isotropic-HytZero gauge effect. At 


R = 250M shown in Fig. 10, the predominant error is the Isotropic-HytZero gauge 
effect. 


These results show that the selected runs do produce a substantial gauge error 
between the worldtube and inertial coordinates and that the Pitt code effectively 












































Figure 4: Waveforms of the real part of the (2, 2) mode of the news function, as computed 
by the Pitt and SpEC codes with extraction worldtube at -R = lOOM. Compared to Fig.[^ at 
this larger extraction radius the worldtube initialization differences lead to a much smaller 
difference between waveforms, which appear nearly identical here. The main discrepancy 
arises from the treatment of the junk radiation at early times. Here, too, gauge differences 
between HytZero and Isotropic are not visible at this scale. 


removes it, while remaining convergent for the duration of the run. 


5.2. SpEC code convergence and removal of gauge effects 

Here we examine the SpEC code’s self convergence for each extraction radii, in the 
same way that the Pitt code was examined in Sec. |5.1[ In Figs. [T^ and [13 


we see 


that convergence, measured with the resolutions indicated in Table is comparable 
to the Pitt code’s convergence, while the potential gauge contamination is consistently 
removed at all worldtube radii. As in Figs. 


and 10, the solid lines (Isotropic) 
and dashed lines (HytZero) overlap due to consistency in gauge removal. The SpEC 
extraction code effectively removes gauge error at all radii while remaining convergent 


















Figure 5: Graph showing the relative difference between the real part of the (2, 2) mode 
SpEC and Pitt news waveforms for extraction radius R = lOOM, and the relative numerical 
error, corresponding to the waveforms in Fig. In comparison with Fig. by i? = lOOM 
the difference between the Pitt and SpEC algorithms has dropped to the level of numerical 
error in each algorithm. 


throughout the runs. 


5.3. Comparison between the Pitt and SpEC codes 


In Secs 0 and |5.2| we have shown that both codes are convergent and remove 
potential gauge effects. We have also demonstrated that the difference in the news 
computed by the two codes disappears as the extraction worldtube radius increases. 
Here we provide further evidence that even at a small worldtube radius the waveform 
computed by the SpEC code is valid. 


After the gauge freedom is removed by extraction, there is still supertranslation 
and Lorentz freedom in the choice of inertial coordinates, which affect the phase and 
velocity of the inertial observers. This effect is highly sensitive to initial conditions 




























































Figure 6: Waveforms of the real part of the (2, 2) mode of the news function, as computed 
by the Pitt and SpEC codes with extraction worldtube at R = 250M. At this large extraction 
radius there is a barely noticeable difference between all the waveforms, limited to the junk 
radiation at early times. 


and also to the evolution of the inertial confomal transformation factor u, especially in 
the extreme gauge conditions of extraction at R = SOM. It feeds into the worldtube 
interpolation error in the Pitt code. In order to verify that the discrepancy illustrated 
in Fig. [^between the news computed by the Pitt null code and SpEC is partially due 
to this inertial coordinate freedom, we compute the time derivative of the news, which 
is related to the Weyl curvature in inertial coordinates according to dtN = This 
suppresses phase differences between the two waveforms. In making the comparisons, 
is computed semi-independently using the Weyl tensor waveform module in the 
current version of the Pitt code [20] • In these runs, was found to be convergent 
with truncation error comparable to the consistency between T® and dtN in the Pitt 
code. 


In Fig. 14, we see that the time derivative of the news and have much less 
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Figure 7: Graphs showing the relative difference between SpEC and Pitt news waveforms at 
R = 250M, in comparison to the relative numerical error, corresponding to the waveforms in 
Fig.H During the post-junk part of the waveform, the error due to worldtube initialization 
has dropped to the level of the numerical error, completing the trend seen in Fig. [5| 


discrepancy than Fig. would suggest. In Figs. 15, 17, and 19, we compare relative 
errors between \[^4 and dtN computed by the Pitt and SpEC codes. Not only is 
there agreement between the codes at R = SOM, this agreement persists for larger 


extraction radii, as shown in Figs. 16 and 18 Both codes show agreement with the 
'F 4 waveform throughout the runs at all three extraction radii. This indicates that a 
major part of the discrepancy in Fig. |^is due to initialization errors in the Pitt code, 
conhrming the physical validity of the extracted SpEC waveform. 


5.4- Relative motion between inertial and worldtube coordinates 

In Section we discussed the construction of an inertial coordinate system and its 
evolution with respect to the worldtube coordinates constructed from the Cartesian 
coordinates of the Cauchy code. Here, we describe the motion of the inertial (0, 0) 




























































































Figure 8: Graphs of the relative error logj^g 1^^22/^221 hr the (2,2) mode of the news 
function for the R = 30M Pitt extraction run in both gauges. The relative errors for 
the Pittl (low) and Pitt2 (high) resolutions (compared to Pitt2 and PittS respectively) 
are rescaled to demonstrate convergence. The dashed blue line indicates the relative error 
(Isotropic vs HytZero) between the news computed in both gauges. The dot-dashed purple 
line (Inertial vs worldtube) indicates the relative error between the news computed in the 
worldtube coordinates and the inertial coordinates. At this small extraction radius, this 
discrepancy is high due to the strong gauge effect of junk radiation. 


angular coordinates relative to the worldtube angular coordinates, constructed in 


the standard way from the worldtube Cartesian coordinates. Figure illustrates 
the global pattern of this relative motion for the Isotropic gauge SpEC run with the 
highest resolution extraction at i? = SOM. Generally speaking, the coordinates wiggle 
back and forth in the direction corresponding to the motion of the black holes. The 
complete movement amounts to at most a few percent of their initial values, but even 
this is sufficient to introduce considerable gauge error in the waveform, as already 
seen in Fig. 
















































Figure 9: Graphs of the relative error log^Q |AA^ 22 /-^ 22 | in the news function for Pitt 
extraction at R = lOOM. Again, the errors for the Pittl and Pitt2 resolutions demonstrate 
convergence. Compared to Fig. the relative error (Inertial vs worldtube) between inertial 
and worldtube coordinates has now dropped below the Isotropic vs HytZero gauge effect. 


In Fig. 21, the relative 0 motion of the point circled in Fig. 20 is plotted as a 
function of inertial time. Initial junk radiation causes considerable wobble, followed 
by a smooth return almost to its starting point. The maximum excursion of the 
^-coordinate shift from its initial value is about 3.5%. 


5.5. Precessing, spinning binary black hole merger 

In addition to the head-on collision tests which we have described, we have also 
investigated stability and convergence of the Pitt and SpEC CCE modules, together 
with the inertial-worldtube coordinate transformation, using the generic test run of 
precessing, spinning binary black holes in [T5], as taken from Taylor et a/. [21]. Its 
parameters are mass ratio g = 3, black hole spins xi = (0.7, 0, 0.7)/a/2 and X 2 = 
(—0.3,0, 0.3)/-\/2, number of orbits 26, total time T = 7509M, initial eccentricity 








































Figure 10: Graphs of the relative error logj^Q | AiV 22 /-/V 22 | in the news function for Pitt 
extraction at i? = 250M. At this large extraction radius, the dominant error arises from 
the Isotropic-HytZero gauge effect. 


10 initial frequency Uini = 0.032/M, and extraction radius R = lOOM. The Pitt 


and SpEC waveforms displayed in Fig. 22 are fairly typical waveforms, spanning the 
initial junk radiation through inspiral, merger, and ringdown. 

Figures and 24 show a log scale comparison of the waveforms with absolute 
error. The codes agree strongly throughout the run. 

The relative inertial-worldtube coordinate motion of a representative point in 


the extended generic run is illustrated in Fig. ^ At early times, the helix has two 
loops per cycle corresponding to each of the black holes. At later times, precession 
dominates the evolution of this particular coordinate. Throughout the run, the 
deviation is around 0. 
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Figure 11: Graphs of the relative error logj^g l^-^ 22 /-^ 22 | in the news function for the 
R = SOM SpEC extraction run in both gauges. The relative errors for the SpECl (low) 
and SpEC2 (high) resolutions (compared to SpEC2 and SpECS respectively) are rescaled 
to demonstrate convergence. The graph Isotropic vs HytZero indicates the relative error 
between the news computed in both gauges. Even at this small extraction radius, there is 
relatively little Isotropic vs HytZero gauge error. 

6. Conclusion 

The SpEC characteristic evolution algorithm has now been furnished with a 
convergent, efficient news extraction module. SpEC is now capable of rapidly 
producing accurate, gauge free waveforms as required. 
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Figure 12: Graphs of the relative error logj^g l^-^ 22 /-^ 22 | in the news function for the 
R = lOOM SpEC extraction run in both gauges. The graphs show convergence in both 
gauges. The Isotropic vs HytZero gauge error is relatively small. 
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Figure 13: Graphs of the relative error logj^Q | AA^ 22 /-^ 22 | in the news function for the 
R = 250M SpEC extraction run in both gauges, showing convergence in both gauges as well 
as small gauge error throughout the run. 
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Figure 16: Waveforms of 'I '4 and dtN as computed for the SpEC and Pitt runs using the 
Isotropic gauge with extraction radius R = lOOM. 
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Figure 17: The relative error log^g between and dtN as computed for 

the SpEC and Pitt runs using the Isotropic gauge with extraction radius R = lOOM, 


corresponding to the waveforms in Fig. 16 Both codes show comparable agreement 
throughout the run. 
































































































0.004 



Figure 18: Waveforms of ^'4 and dtN as computed for the SpEC and Pitt codes using the 
Isotropic gauge with extraction radius R = 250M. 



Figure 19: The relative error log^g between T 4 and dtN computed for 

the SpEC and Pitt runs using the Isotropic gauge with extraction radius R = 250M, 


corresponding to the waveforms in Fig. 18 Both codes show good agreement throughout 


the run. 
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Figure 20: The arrows indicate the global pattern of the motion of the (0, 0) inertial 
coordinates at each spherical collocation point (labeled in radians) relative to the worldtube 
angular coordinates. Global motion is primarily in the Cartesian x direction, and has been 
exaggerated by a factor of 10 for visibility. 
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Figure 21: The relative inertial-worldtube (/)-coordinate motion of the point circled in 


Fig. 20 shows approximately a 3.5 percent variation from its initial value. 
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Figure 22: Waveform of the real part of the (2, 2) mode of the news function for the generic 
processing binary black hole run, showing 26 orbits, inspiral, merger, and ringdown. 
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Figure 23: Graphs of log^Q |i?e(news 22 )| for Pitt and SpEC waveforms for the processing 
binary black hole run, showing 16 orbits, inspiral, merger, and ringdown. The difference 
Pitt — SpEC gives the absolute error log;^o \^22piu ~ ^22spEc\ between the codes, showing 
consistent agreement throughout the run. 








































































































































Figure 24: 


Same graphs as in Fig. 
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but focusing on the merger part of the waveform. 


Absolute error shows consistency between the Pitt and SpEC news waveforms through 


merger and ringdown. 
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Figure 25: Relative motion of a representative angular coordinate taken from SpEC3 with 
extraction radius R = lOOM showing a long term helical pattern concordant with the 
black hole inspiral and merger. Globally, initial oscillation due to junk radiation is aligned 
primarily along the Cartesian x direction, corresponding to the initial black hole orientation. 
Coordinate motion is an epicyclic helix whose amplitude is modulated by precession of the 
orbital plane. 
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